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INTRODUCTION 


Typical  electromagnetic  shielding  problems  involve  the  calcu- 
lation of  the  time-dependent  electric  or  magnetic  field  penetration  of 
electrically  conducting  media.  Two  material  properties  which  enter 
into  such  calculations  are  the  electrical  conductivity  0 and  the 
permeability  y.  In  most  cases  the  electrical  conductivity  does  not 
significantly  vary  with  electromagnetic  field  level  and  can  be  re- 
garded as  a constant  for  a given  material.  In  nonmagnetic  materials 
the  permeability  is  also  independent  of  field  level  and  has  the  value 
for  free  space  (vacuum)  Po  = Att  x 10“^H/m.  For  media  whose  material 
properties  are  independent  of  the  applied  electromagnetic  field  level, 
the  partial  differential  equations  which  govern  electromagnetic  field 
propagation  are  linear.  The  analytical  theory  for  electrc^nagnetic 
field  calculations  in  media  with  constant  permeability  has  been  well 
developed  and  extensive  analyses  exist.  In  fact,  quite  often  the 
exact  analytical  solution  can  be  achieved  by  any  one  of  a number  of 
standard  linear  mathematical  techniques.  ' 

On  the  other  hand,  ferromagnetic  materials,  which  are  common- 
ly used  in  electromagnetic  shield  fabrication,  not  only  have  a permea- 
\ bility  that  is  dependent  on  the  magnetic  field  strength  H but  also 

I exhibit  the  phenomenon  of  saturation.  Although  ferromagnetic 

: materials  usually  have  relatively  large  (in  some  cases  relatively 

■ enormous)  permeabilities  compared  to  nonmagnetic  materials,  the  perme- 

I ability  can  be  greatly  reduced  if  the  material  is  driven  into  satura- 

i tion.  For  these  materials  the  partial  differential  equations  which 

i govern  electromagnetic  field  propagation  are  nonlinear  and 


CROISANT  & NIELSEN 


ET 


considerably  less  amenable  to  analytic  solution.  Due  to  the  difficul- 
ties associated  with  the  nonlinear  partial  differential  equations 
which  arise  in  this  case,  relatively  few  analyses  are  available  so 
that  it  is  difficult  to  assess  the  effect  of  a variable  permeability 
on  the  time-dependent  electromagnetic  field  penetration. 

This  study  examines  the  time-dependent  pepetration  of  a step 
increase  in  magnetic  field  strength  H into  a semi-infinite  conducting 
medium  having  a permeability  which  varies  with  magnetic  field 
strength.  The  medium  is  considered  to  be  isotropic,  homogeneous,  and 
Initially  demagnetized.  The  medium  at  any  point  is  presumed  to  follow 
its  initial  magnetization  curve  for  which  a simple  approximation  is 
assumed.  The  analytical  approach  consists  of  mathematical  analysis 
supplemented  with  numerical  calculations.  The  problem  can  be  simpli- 
fied from  one  involving  a second  order  nonlinear  partial  differential 
equation  to  one  involving  a second  order  nonlinear  ordinary  differen- 
tial equation  utilizing  a simple  transformation  of  variables.  Al- 
though a formal  parametric  representation  of  the  solution  of  the 
second  order  differential  equation  is  considered,  a simple  closed  form 
solution  in  terms  of  elementary  functions  does  not  seem  to  exist. 

Hence  several  analytical  techniques,  as  well  as  numerical  calcula- 
tions, have  been  employed  to  deduce  properties  of  the  solution. 

MATHEMATICAL  FORMULATION  OF  THE  PROBLEM 


The  propagation  of  electromagnetic  fields  in  a conducting 
medium  (where  the  displacement  current  can  be  neglected)  is  governed 
by  the  equations 


V X E = - 


at 


(1) 


and 


V X H = 0E  , 


(2) 


where  E is  the  electric  field,  H is  the  magnetic  field  strength,  and 
S is  the  magnetic  flux  density.^  In  a_^homogeneous  and  isotropic 
medium,  it  can  be  assumed  that  B = B(H) , i.e.,  that  the  magnetic  flux 
density  is  a function  of  the  applied  magnetic  field  strength  and  is  in 
the  same  direction.  With  the  differential  permeability  defined  as 
IJd(H)  = dB/dH,  (1)  can  be  written  as: 

V X E - - . (3) 

Consider  a homogeneous,  isotropic,  conducting  medium  of  semi- 
infinite  extent  (x  ^ 0) . In  the  present  problem  it  is  presumed  that 
the  medium  is  initially  demagnetized,  and  at  time  t = 0,  a step  in- 
crease in  magnetic  field  strength  in  the  y direction  occurs  at  the 
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surface  x = 0,  e.g. , as  the  result  of  an  imposed  surface  current  in 
the  z direction. 


Since  it  is  presumed  that  H = Hy(x,t)y,  the  propagation  of 
the  electromagnetic  fields  following  the  onset  of  the  step  increase  in 
Hy(x,t)  is  governed  by  the  equations 


8E^(x,t) 


8t 


3H  (x,t) 

-L—  • aE^Cx.t)  , (5) 

By  eliminating  E2(x,t)  from  (4)  and  (5),  it  is  found  that  the  partial 
differential  equation  governing  the  penetration  of  H^Cx.t)  is 

9^H  (x,t)  3H  (x,t) 

-fr-  = • <^> 

9x 

In  general,  B(H)  = where  M(H)  is  the  magnetization.  As 

H is  increased  from  zero,  usually  starts  from  an  initial  value 

Pi  is  usually  larger  than  Uq) , increases  to  some  maximum  value, 
and  then  decreases  to  yg  as  the  material  undergoes  saturation.  In 
many  cases  before  saturation,  the  variation  of  B with  H [hence  the 
variation  of  yj}(H)]  is  primarily  due  to  M(H) . In  this  study  it  is 
assumed  that  over  the  range  of  interest  B(H)  can  be  represented  by  the 
simple  approximation 

B(H)  = Bg[l  - exp(-H/Hj], 

where  Bg  represents  a saturation  value  for  B and  is  a parameter 
which  Indicates  the  shape  of  the  magnetization  curve  for  the  medium. 
With  this  representation  of  the  magnetization  curve, 

y^(H)  = y^  exp(-H/H^)  , 

where  y^^  = Bg/H^  is  the  initial  slope  of  the  B-H  curve.  The  variation 
of  B(H)  with  H is  shown  in  Figure  la  and  the  corresponding  variation 
of  yd(H)  with  H is  shown  in  Figure  lb.  For  large  H,  y^jCH)  actually 
approaches  yg  as  a lower  limit;  hence  this  simple  representation  is 
appropriate  as  long  as  y^^  exp(-H/Hjj)  > y^  over  the  range  of  H under 
consideration. 
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FIGURE  la  B(H)  VS.  H 


H/Hm 

FIGURE  lb  vs.  H 


The  present  problem,  therefore,  requires  the  solution  of  the 
nonlinear  partial  differential  equation 


3 H (x,t) 
. 2 


= ay^  exp(-H/H^) 


3H  (x,t) 
y 

3t 


subject  to  the  initial  condition 

Hy(x,0)  = 0 

and  the  auxiliary  conditions 

H (0,t)  = H 
y o 


HyC^.t)  = 0 . (10) 

The  electric  field  can  be  determined  from  the  solution  to  the  above 
problem  by  noting  from  (5)  that 

, 3H  (x,t) 

= a 3x  " • 

REDUCTION  IN  THE  NUMBER  OF  INDEPENDENT  VARIABLES 


Using  an  approach  similar  to  that  used  for  nonlinear  diffu- 
sion phenomena  analysis,  the  number  of  independent  variables  appearing 
in  the  problem  can  be  reduced  by  the  transformation 
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Hy(x,t)  = H(C); 


? 


i x_ 
2 /t 


The  new  variables  are  introduced  by  evaluating  the  partial 
derivatives  in  terms  of  the  new  variables  and  substituting  these 
quantities  into  the  original  partial  differential  equation.  On 
calculating  the  partial  derivatives  it  is  found  that 

3H^(x,t)  ^ ^ dH(^)  ^ dH(C)  , 

3x  3x  dC  2/t” 


and 


3‘‘h  (x.t) 
3x^ 


^ ^ dH(^)  /35^\" 

5x2  d?  Vax; 


2 d^H  (O 


*^1  d^H(^) 


4t 


dr 


3Hy(x»t)  ^ ^ ^ _ 1 


'o\i. 


3t 


3t  di; 


2 2 ,3/2  dC 


dH(C)  ^ _ 1 1 dH(i;) 


2 t dC 


Upon  substituting  these  quantities  into  (7)  and  cancelling  a factor  of 
ay^/4t,  it  is  found  that  the  magnetic  field  strength  satisfies  the 
equation 

“ = - 2C  exp(-H/H  ) , (12) 

d^2  m d? 

in  which  the  variables  x and  t no  longer  appear  explicitly.  The 
electric  field  is  related  to  H(?)  by 


E^(x,t)  = 2 


i dH(^) 

a ^ dC 


(13) 


in  which  a factor  of  ^ appears.  Noting  that  both  x = 0 and  t = “ 
imply  that  ^ = 0 and  that  both  x = <»  and  t = 0 imply  that  C = it 
follows  from  (9)  that 

H(0)  = H (14) 

o 

and  from  (8)  and  (10)  that 

Hy(<»)  = 0 . (15) 

For  computational  purposes,  it  is  convenient  to  normalize  the 
problem  in  terms  of  H . With 
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F(C)  = H(C)/H^  = Hy(x,t)/H^  , 

(16) 

the  problem  has  been  reduced  to  solving 
ordinary  differential  equation 

the  nonlinear  second  order 

= - 2^  exp(-aF)  ^ 

AX,  ^ 

(17) 

subject  to  the  conditions 

F(0)  = 1 

(18) 

and 

F(<»)  = 0 , 

(19) 

where  a = Hq/H^.  Having  solved  for  F(^) , the  electric  field  can  then 
be  evaluated  from  the  relation 


dF(;:)  _ ^ /o_ 

« ■ “,  / 


/T  E (x,t)  . 
z 


ANALYTICAL  ASPECTS  OF  THE  ORDINARY  DIFFERENTIAL  EQUATION 


(20) 


Classical  Linear  Solution 


In  classical  analyses  of  electromagnetic  field  penetration, 
it  is  assumed  that  the  permeability  is  constant,  which  corresponds  to 
a = 0 in  the  present  investigation.  For  a = 0 the  problem  reduces  to 


(21) 


subject  to  conditions  (18)  and  (19).  This  equation  can  be  integrated 
in  a straightforward  manner,  to  obtain  the  classic  solution 

F(C)  = 1 - erf(C)  , (22) 


where  the  error  function  is  defined  as 


erf(i;) 


— / exp(-C^)dCj^  . 

/n  0 


It  follows  Immediately  that 

“ exp(-?^) . 

/iT 


dF(;) 

dC 


(23) 


I 


» 
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Formal  Solution  of  the  Nonlinear  Problem 


Most  of  the  standard  linear  techniques  are  of  relatively  lit- 
tle use  in  achieving  an  exact  solution  in  nonlinear  cases  (a  0) ; 
however,  in  the  course  of  an  investigation  of  a related  problem  in  non- 
linear diffusion  phenomena,  Fujita  [1]  found  a formal  solution  which 
can  be  directly  adapted  to  the  problem  presently  under  investigation. 
The  problem  is  satisfied  with  F and  ^ being  given  by  the  parametric 
representations 


exp 


3Jin(C  ) 


1/2 

-1/2 


- 3£n(C^)^  d^A  , 


(24) 


(25) 


where  the  derived  quantity  B is  related  to  the  given  quantity  a by 


a = 2 


- ei^n(?J) 


(26) 


The  parametric  representation  of  dF/d4  can  be  evaluated  from 


iE 

d4 


- 


(27) 


The  parameter  C (0  <.  C i D relates  a value  of  F to  the 
corresponding  value  of  (Note  that  s = 0 implies  ;;  = <»  and  F = 0, 
while  5=1  implies  4=0  and  F = 1.)  To  evaluate  F versus  4,  a value 
of  5 Is  selected  and  the  value  of  F and  the  corresponding  value  of  4 
are  evaluated.  The  integrals  appearing  in  the  formal  solution  do  not 
appear  to  be  evaluated  in  closed  form,  thus  numerical  integration 
seems  to  be  required;  howevei;  the  above  representation  offers  the 
advantage  of  allowing  the  computational  effort  to  be  directed  to  the 
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very  accurate  calculation  of  the  solution  at  only  a few  points,  or  one 
point  for  that  matter,  since  the  knowle.^^-  of  intermedi.  ze  points  is 
not  required.  Of  special  interest  is  dF/dC  evaluated  at  C = 0 where 

^ exp(-a/2)  . (28) 


dF(0)  _ 
dC 


aB 


1/2 


I.  ddition  to  its  relation  to  the  electric  field  at  x = 0,  the  cal- 
iated  value  of  dF(0)/dC  can  be  used  to  transform  the  problem  from  a 
two-point  problem  (with  a condition  at  x = “ as  well  as  one  at  x = 0) 
to  an  initial  value  problem  (with  both  boundary  conditions  specified 
at  X = 0) . Such  a transformation  is  often  useful  from  a computational 
standpoint. 

Taylor  Series  Expansion  of  the  Solution 

If  F(^)  was  a known  function  of  ^ and  if  F(^)  and  its 
derivatives  were  continuous,  then  F(C)  could  be  expanded  in  a Taylor 
(Maclaurin)  series  about  ? = 0 


HO  = I 

m=0 


d"*F(0) 
dc“ 


L . F(o,  * iSM 


iC, 


where  F and  its  derivatives  are  evaluated  C = 0.  For  example,  the 
solution  (22)  to  the  linear  case  (a  = O'  has  the  series  representation 


ViO  = 1 - 


I 


^n=0  <2n+l)nl 


2 


21 


- ...) 


A series  expansion  approach  can  also  be  used  to  develop  the 
solution  F(C)  from  the  second  order  differential  equation.  If  F(0) 
and  dF(0)/dC  are  given  as  boundary  conditions,  then  the  approach  is 
straight-forward  since  the  second  and  higher  order  derivatives  can,  in 
principle,  be  determined  from  successive  differentiation  of  the  second 
order  ordinary  differential  equation.  Although  F(0)  is  one  of  the 
auxiliary  conditions  in  the  problem  presently  under  investigation,  it 
is  evident  that  a difficulty  arises  because  the  second  condition  F(») 

= 0 is  specified  at  C = “•  For  this  reason  we  are  led  to  consider 
the  associated  initial  value  problem 


= - 2C  exp(-aF)  ^ 


subject  to 

F(0) 


(29) 


(30) 
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and 


dF(0) 

d? 


= Y • 


(31) 


where  y is  to  be  chosen  such  that  F(“)  = 0.  This  value  for  y can  be 
computed  from  Eq.  (27). 

The  second  derivative  can  be  determined  directly  from  (29) 
from  which  it  follows  that 


d^F(O) 

dC^ 


= 0 . 


The  third  derivative  can  be  found  by  differentiation  of  Eq.  (29) 
^ = - 2 |exp(-aF)||  + [exp(-aF)  g]  | , 


(32) 


from  which  it  follows  that 


d F(0)  ^ 
dC^ 


2 y exp (-a)  . 


Differentiating  (32)  and  evaluating  the  result  at  ^ = 0 yields 

d^F(O)  „ 2 , . 

7 — = -2ay  exp (-a)  . 

dr 

Similarly,  it  is  found  that 

= -6a^y^  exp  (-a)  + 12y  exp  (-2a)  . 

dr 

The  approach  could,  in  principle,  be  extended  Indefinitely;  however, 
in  practice  it  becomes  rather  tedious  after  the  first  few  terms. 

To  fifth  order  terms  the  solution  F(^)  is  given  by 


F(r  = 1 + yc  - 


2 3 


, 2y  exp (-2a)  -ay  exp(-a)  ,5 
20  ^ 


(33) 
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The  series  representation  (33)  can  be  used  to  calculate  F(C)  for  small 
values  of  C. 

In  the  study  of  pulse  penetration  it  is  frequently  of 
interest  to  know  the  time  that  it  takes  for  the  fields  at  a certain 
point  to  rise  to  a given  value.  For  this  purpose  a convenient 
manipulation  is  the  reversion  of  series.  Given  a series  for  the 
dependent  variable  y in  terms  of  the  independent  variable  x, 

y = ax  + bx  + cx  + dx  + ex  + . . . . 

Using  a reversion  of  series  one  can  write  x as  a function  of  y: 


where 


x = Ay  + By^  + Cy^  + Dy^  + Ey^  + . . . , 


b 2b  -ac  „ 5abc  - a d - 5b 

Aa»°  3’^”  5 7 

a a a 

„ 6a^bd  + 3a^c^  + 14b^  - 21ab^c  - a^e 


For  example,  the  classic  linear  solution  (22)  can  be  rewritten  as 


= 1.  y (-1)“  .2n+l  _2_  c!  . ci  ^ 

A reversion  of  series  yields 

C - [ (1  - F)  + (I  - F)3  + 11^  (1  - F)5  + ...  ] 

which  is  suitable  for  calculating  the  value  of  C at  which  the  solu- 
tion F(C)  has  a specified  value  for  small  values  of  (1  - F)  (which 
implies  small  ?). 

Equation  (33)  can  be  rewritten  as 

(1  - F)  . -YC  -F  X-Etal  ,.3  . 2^  e«p(-a)c'' 


2y  exp (-2a)  -ay  exp (-a)  5 

20  C + ...  . 
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yields 

F)  + _ f)3  + CLexE^zoO  ^ _ ^^4 

3y  6y 

A a%xp(-a)  7 exp  (-2a)  , , .5 

'•20  2 30  4 ; (.X  r;  -r  . . . 

Y Y 

which  can  be  used  [for  small  values  of  (1  - F) ] to  examine  the  varia- 
tion with  a of  the  value  of  C at  which  F(^)  has  a specified  value. 

NUMERICAL  CALCULATIONS 

For  a specified  value  of  B the  corresponding  value  of  a can 
be  computed  by  numerical  integration  of  (26)  and-  the  value  of  dF(0)/dC 
= Y can  be  calculated  from  (28).  Some  numerical  calculations  were 
performed  and  the  results  are  shown  in  Table  1.  An  examination  of 
Table  1 reveals  that  rather  large  changes  in  6 correspond  to  relative- 
ly small  changes  in  a.  This  situation  can  make  it  difficult  to 
determine  the  value  of  3 which  corresponds  to  a specified  value  of  a. 
On  the  other  hand,  it  is  evident  that  dF(0)/dC  = Y exhibits  a 
relatively  moderate  variation  with  a.  l-fhile  a varies  in  the  range 
0 s a < 0°  , Y varies  in  the  range  - 2//n  < y S 0.  This  moderate 
variation  makes  it  possible  to  accurately  represent  Y as  a function 
of  a using  the  regression  equation 

- 1.1284  + 2.0490  x lO'^a  - 3.4426  x 10"^a^ 

+ 4.4139  X 10" V - 4.1562  x 10"‘^a'^ 

+ 2.7025  X 10"\^  - 1.0736  x 10"^a^ 

+ 0.9395  X 10"^a^  - , 

which  provides  a convenient  means  of  evaluating  y for  0 < a < 10. 

The  variations  of  F and  dF/d?  with  Z were  calculated  for  a = 
0(3  = °“) , 01  = 2(3  = 3.959  x 10“^),  and  a = 5(3  = 5.922  x 10“3)  and  are 
shown  in  Figures  2a  and  2b.  The  magnetic  field  strength  H can  be 
determined  from  (16)  and  the  electric  field  from  (20) . In  shielding 
applications,  the  time  variation  is  often  of  interest;  therefore  the 
variations  of  F and  dF/dC  with  1/^2  are  shown  in  Figures  3a  and  3b. 


dF(Q)  ^ 
dC 


A reversion  of  series 

5 u r (1  - 
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Table  1 

Numerical  Calculations  of  a and  dF(0)/dC.* 


dF(0) 

dC 


00 

0 

-Z//n 

5 X 10“^ 

0.0112 

-1.1260 

1 X 10^ 

0.0249 

-1.1233 

5 X 10^ 

0.0351 

-1.1212 

1 X 10^ 

0.0773 

-1.1127 

5 X 10^ 

0.1083 

-1.1066 

1 X 10^ 

0.2326 

-1.0825 

5 X 10^ 

0.3197 

-1.0662 

1 X 10^ 

0.6424 

-1.0098 

5 X 10*^ 

0.8475 

-0.9769 

1 X 10° 

1.5058 

-0.8847 

5 X 10“^ 

1.8681 

-0,8414 

1 X lO”^ 

2.8683 

-0,7432 

5 X lO”^ 

3.3552 

-0.7043 

1 X 10“^ 

4.5797 

-0.6255 

5 X lO"^ 

5.1378 

-0.5965 

1 X 10"° 

6.4838 

-0.5392 

—A 

5 X 10  ^ 

7.0801 

-0.5183 

-4 

1 X 10 

8.4931 

-0.4767 

5 X 10"^ 

9.1115 

-0.4612 

1 X 10"^ 

10.5651 

-0.4300 

*The  calculations  were  performed  on  a CDC  6600  computer  using  simple 
trapezoidal  numerical  integration. 


I 
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CONCLUSIONS 

In  the  linear  case  (a  = 0),  the  solution  F(?)  is  linearly 
scaled  by  Hq  and  the  shape  of  the  solution  F(C)  remains  unchanged. 
Examination  of  the  preceding  analyses  and  numerical  calculations 
reveals  that  the  nonlinear  case  is  more  complicated.  As  might  be 
expected,  the  results  in  the  nonlinear  case  may  not  be  simply  scaled 
in  the  usual  linear  fashion.  Although  by  the  nature  of  the  problem 
the  final  value  Hq  is  the  same,  the  manner  in  which  this  value  is 
approached  depends  on  the  value  of  Hq.  The  variation  in  permeability 
with  applied  magnetic  field  strength  in  the  present  problem  results 
in  a more  rapid  penetration  of  the  medium  by  the  magnetic  field.  In 
other  words,  at  a given  position  in  the  medium,  the  magnetic  field 
strength  H will  reach  a given  percentage  of  the  applied  field  Hq  in  a 
shorter  time.  Alternatively,  at  a given  time  the  point  at  which  the 
magnetic  field  strength  has  a prescribed  value  will  have  penetrated  to 
a greater  distance.  In  addition  it  is  evident  that  for  a given  value 
of  the  parameter  a the  solution  F(C)  at  some  ^ = c will  have  the  same 
value  for  all  combinations  of  x and  t which  are  related  by  i;  = c. 

This  indicates  that  if  a certain  value  for  the  magnetic  field  strength 
H has  reached  location  x at  some  time  t,  then  the  same  value  of  H will 
occur  at  some  location  x'  = ax  at  some  time  t'  = a^t.  It  might  be 
pointed  out  that  for  a given  a,  the  solution  F(^)  will  have  reached  a 
prescribed  value  Fp  at  some  ^p.  Thus,  one  can  consider  a "penetration 
thickness" 

6 5 c (2/^ir)/r 
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beyond  which  the  magnetic  field  strength  has  changed  by  less  than  'p' 
percent  of  H^.  For  example,  the  linear  solution  (22)  has  a value  less 
than  0.01  when  the  error  function  argument  is  about  2.  If  it  is 
necessary  to  calculate  the  electromagnetic  fields  in  a shield  of 
finite  thickness,  then  the  solution  F(C)  will  be  a good  approximation 
if  the  "penetration  thickness"  is  small  with  respect  to  the  shield 
thickness. 
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